Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

Role of Equatorial Anomaly in Earthquake time precursive features: A few strong events over West Pacific zone

The earthquake (EQ) time coupling processes between equator-low-mid latitude ionosphere are complex due to inherent dynamical status of each latitudinal zone and qualified geomagnetic roles working in the system. In an attempt to identify such process, the paper presents temporal and latitudinal variations of ionization density (foF₂) covering 45°N to 35°S, during a number of earthquake events (M?>?5.5). The approaches adopted for extraction of features by the earthquake induced preparatory processes are discussed in the paper through identification of parameters like the ‘EQ time modification in density gradient’ defined by δ?=?(foF_{2 max}???foF_{2 min})∕τ_mm, where τ_mm – time span (in days) between EQ modified density maximum and minimum, and the Earthquake time Equatorial Anomaly, i.e. EEA, one of the most significant phenomenon which develops even during night time irrespective of epicenter position. Based on the observations, the paper presents the seismic time coupling dynamics through anomaly like manifestations between equator, low and mid latitude ionosphere bringing in the global Total Electron Content (TEC) features as supporting indices.     

Diurnal and seasonal variations in long-term changes in the E-layer critical frequency

Long-term changes in the E-layer critical frequency, foE, at three stations of the European region (Juliusruh, Slough and Rome) and also at Moscow and Wakkanai stations are analyzed by the method developed by the authors and described in detail in the previous papers. It is found that Juliusruh and Slough stations demonstrate a well-pronounced change in foE (a trend) during two previous decades. At the same time, the same features of the behavior of the aforementioned trend k(foE) are obtained. The trend is positive and negative in the morning and evening hours, respectively. Similar diurnal behavior of k(foE) is found also for Moscow station but with lower absolute values of the trends. A well-pronounced seasonal behavior of k(foE) is detected at Juliusruh and Slough: the trend is minimal and maximal in the summer period and at the end of fall—beginning of winter, respectively. The maximal amplitude in the morning hours reaches +0.04?MHz per year, whereas the minimal amplitude in evening hours is ?0.06?MHz per year. No systematic changes exceeding by the magnitude 0.01?MHz per year are found for Rome and Wakkanai stations. It is assumed that the observed trends are related to changes (trends) in the meridional wind bringing NO molecules from the auroral oval to lower latitudes.     

Assessment study of ionosphere correction model using single- and multi-shell algorithms approach over sub-Saharan African region

It is a known fact that ionosphere is the largest and the least predictable among the sources of error limiting the reliability and accuracy of Global Navigation Satellite Systems (GNSS) and its regional augmentation systems like Satellite Based Augmentation System (SBAS) in a safety-of-life application. The situation becomes worse in the Equatorial Ionization Anomaly (EIA) region, where the daytime ionization distribution is modified by the fountain effect that develops a crest of electron density at around ±15° to ±20° of the magnetic equator and a trough at the magnetic equator during the local noon hours. Related to this phenomenon is the appearance of ionosphere irregularities and plasma bubbles after local sunset. These may degrade further the quality of service obtained from the GNSS/SBAS system of the said periods. Considering the present operational augmentation systems, the accuracy and integrity of the ionosphere corrections estimate decreases as the level of disturbances increases. In order to provide a correct ionosphere correction to the user of GNSS operating in African EIA region and meet the integrity requirements, a certified ionosphere correction model that accurately characterizes EIA gradient with the full capacity to over-bound the residual error will be needed. An irregularities detector and a decorrelation adaptor are essential in an algorithm usable for African sub-Saharan SBAS operation. The algorithm should be able to cater to the equatorial plasma vertical drifts, diurnal and seasonal variability of the ionosphere electron density and also should take into account the large spatial and temporal gradients in the region. This study presents the assessment of the ionosphere threat model with single and multi-layer algorithm, using modified planar fit and Kriging approaches.     

Precise line-of-sight modelling for angles-only relative navigation

This work presents a precise analytical model to reconstruct the line-of-sight vector to a target satellite over time, as required by angles-only relative navigation systems for application to rendezvous missions. The model includes the effects of the geopotential, featuring: the analytical propagation in the mean relative orbital elements (up to second-order expansion), the analytical two-way osculating/mean orbital elements’ conversion (second-order in $J_2$ and up to a given degree and order of the geopotential), and a second-order mapping from the perturbed osculating elements’ set to the local orbital frame. Performances are assessed against the line-of-sight reconstructed out of the precise GPS-based positioning products of the PRISMA mission. The line-of-sight modelled over a far-range one day long scenario can be fitted against the true one presenting residuals of the order of ten arc-seconds, which is below the typical sensor noise at far-range.     

Atmospheric reentry hemisphere prediction for prograde orbits using logical disjunction

A unique logic-based algorithm for atmospheric reentry hemisphere prediction is presented for spacecraft in low-eccentricity, prograde low Earth orbits at altitudes of 300 km and lower. Using two-line element (TLE) data for initial orbit conditions, coupled with coarse estimates for spacecraft aerodynamic characteristics, the algorithm relies on logical disjunction operations based on a dual analysis of histogram and two-weighted Gaussian probability density function (PDF) fits of predicted reentry latitude data. The algorithm requires the execution of a series of parametric simulations to determine the reentry hemisphere for variations in spacecraft aerodynamic coefficients and drag reference area. When implemented, the algorithm yields accurate hemisphere predictions on average 15 days from reentry as demonstrated by historical reentry cases from 1979 to 2018. All reentry cases were selected to demonstrate the algorithm’s ability to deliver accurate reentry hemisphere predictions for spacecraft with varying physical size and mass, and reentering during different periods of solar cycle activity.     

Doppler characteristics of HF spread Doppler clutter recorded at Alice Springs,Australia

Mean night-time peak power, Doppler shift and Doppler width of spread Doppler clutter (SDC) received by a high frequency backscatter radar located at Alice Springs, Australia from 2000 to 2018 is presented as a function of azimuth, sunspot number, time of year and frequency. The sampled region covers 90 degrees from West to North and includes the northern and southern equatorial anomalies.SDC peak power diminished across all azimuths during the winter solstice from around May to August (local winter) coinciding with the global decrease in F layer density due to the annual non-seasonal F2 anomaly but was generally constant during the equinoxes. In contrast, SDC Doppler width and inbound Doppler shift both increased during the equinoxes and exhibited azimuthal dependence related to the eastward equatorial plasma drift.SDC peak power increased with increasing sunspot number with frequency dependence during winter but not summer. Inbound Doppler shift and Doppler width increased with increasing sunspot number during equinox but not solstice with a strong dependence on azimuth and a weak dependence on frequency.     

Study of highly perturbed spacecraft formation dynamics via approximation

This paper explores methods for approximating and analyzing the dynamics of highly perturbed spacecraft formations with an emphasis on computationally efficient approaches. This facilitates on-board computation or rapid preliminary mission design analysis. Perturbed formation dynamics are often approximated as linear time-varying (LTV) systems, for which Floquet theory can be used to analyze the degree of system instability. Furthermore, the angular momentum of the relative orbital state can be computed with the approximate dynamics to provide additional insight. A general methodology is developed first and then applied to the problem of unstable formation dynamics in asteroid orbits. Here the dominant perturbative effects due to low-order gravitational harmonics and solar radiation pressure are modeled. Numerical simulations validate the approach and illustrate the approximation accuracy achieved.